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Abstract 

We study the effects of weak long-ranged antiferromagnetic interactions of 
strength Q on a spin model with predominant short-ranged ferromagnetic 
interactions. In three dimensions, this model exhibits an avoided critical 
point in the sense that the critical temperature T C (Q = 0) is strictly greater 
than limQ^oT c (Q). The behavior of this system at temperatures less than 
Tc(Q = 0) is controlled by the proximity to the avoided critical point. We 
also quantize the model in a novel way to study the interplay between charge- 
density wave and superconducting order. 
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I. INTRODUCTION 



There are many diverse contexts in which a short-ranged tendency to order is thwarted 
by a long-ranged frustrating interaction. In particular, recent theories of the glass transition 
i l[ and of the properties of doped Mott insulators [|2],|3| have lead to the consideration of 
such models. In the first example, a tendency of the molecules in a supercooled liquid to 
pack into a locally preferred structure is frustrated by the inability of such structures to tile 
space; the long-range nature of the induced interaction represents the superextensive growth 
of strain which would occur if the locally preferred structure were forced to tile space. In the 
second example, a short-range tendency of holes in an antiferromagnet to phase separation 
competes with the long-range Coulomb repulsion between holes. In this latter case, effects 
of quenched disorder and quantum fluctuations may also be important. 

It has been argued fl[ that the effect of such uniform, long-range frustration is an avoided 
critical point, leading to a phase diagram of the sort shown in Figure 1. Here, T is the 
temperature and Q is the strength of the frustrating interaction. The salient feature of this 
phase diagram is the discontinuity in T c in the Q — > limit. The presence of the avoided 
critical point leads to new types of fluctuation phenomena for Q ^ 1 and T < T C (Q = 0). 

To capture the essential physics of these problems, we will consider a spin model on the 
(i-dimensional hypercubic lattice. Our starting point is the classical Hamiltonian 

H cl [S}= J - £ \S n - ^,| 2 + | E v(R - R')S n ■ Sg, (1.1) 

<R,R'> R¥=R' 

where J and Q are positive couplings, the sum over < R,R' > is over nearest neighbor pairs, 
and v(R) is long-ranged, with asymptotic falloff 

v(R) ~ \R\~ X (1.2) 

with < x < d. Specifically, it is most convenient to express the interaction on the lattice 
in terms of its Fourier transform, which in the simplest cases we will consider here, implies 
that the Hamiltonian is 
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H c i[S] = Wj$) V- (1-3) 
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Here 5* is the Fourier transform of S, 



J(k) = JA(k)+Q[A(k)] y + ... (1.4) 

with y = (d — x) /2, and 

A$=£;V(*a) (1.5) 
a=l 

is the Fourier transform of the lattice Laplacian, and 

V(k) = 2 [1 - cos(Jfe)] « A; 2 for jfe < 1. (1.6) 

In Eq. fln| ), the ... refers to additional, short-range terms which are generically present and 
which will be included for reasons that will become clear shortly. To make our discussion 
explicit, and because it is typically the case of most physical interest, we will focus on the 
Coulombic case y = 1; the principal results are qualitatively similar for any long-range 
potential. 

This model can be studied with various definitions of the spin variables, Sg. Probably, 
the most interesting case is the Ising version, [§-§1 where Sg = ±1. The majority of the 
results in the present paper concern the exactly solvable "mean spherical" version of this 
model, in which the spins are taken to be real numbers with the mean global constraint 

N- ly £<[S^ 2 >=l (1.7) 

R 

where N is the number of sites. As is well known, this model is equivalent, in the ther- 
modynamic limit, to the usual spherical model in which this global constraint is enforced 
configuration by configuration J7| and to the large n limit of the 0(ri) model. |J (Ongoing 
work on the 1/n expansion will shed some additional light on the relation of the present 
results to the properties of the model at finite n. 0) 

Since, in the context of doped Mott insulators, quantum effects are important, we also 
consider a quantum version of this model. In particular, we wish to consider a model with 



two distinct types of low temperature ordered phases: those with spin (charge) order and 
those with momentum (superconducting) order. These order parameters are canonically 
conjugate and dual to each other. We are motivated in this choice by an analogy with 
models of hard-core bosons on a lattice, or equivalently to a spin s = 1/2 frustrated quantum 
Heisenberg-Ising model in which the XY coupling is the particle kinetic energy (and hence 
associated with the momentum) while the Ising ordering is associated with ordering of the 
dual fields, (charge ordering) which can either correspond to phase separation or formation 



of a charge density wave ordered state. [LOj With this in mind, we define a quantum model 



H = H cl [S)+H qu [P], (1.8) 

where H qu and H c i are appropriate quadratic forms, in which the "momenta", {-P^}, are 
cannonically conjugate to the spins, {S^;}, i.e. 

[S n ,P n ,]=ih5 M)n ,. (1.9) 

The constraint equation is 

N- 1 Y.[*s(\S n \ 2 ) + a p (\P A \ 2 )}=s 2 . (1.10) 

R 

So long as neither a s nor a v is zero, we can, without loss of generality, [[□]] take ot s = a p = 1. 



To be concrete, we also confine our considerations to the simplest case in which H qu consists 
of the simplest, unfrustrated nearest-neighbor ferromagnetic interaction 

H qu [P] = Y £ \ P A- P Ri\ 2 - (I-") 

<K,-R'> 

Thus, W is the bare superfluid stiffness. It also is necessary to augment the Hamiltonian 
with a uniform field, K, which favors charge ordering (< Sg >^ 0) when negative and 
superconducting ordering ((< Pg >^ 0) when positive, 

H c1 ^H c1 + ^Y,\ S r\^ (1-12) 
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R 



or, in other words, J{k) — > J{k) + K. 



In the context of quantum spin glasses, Nieuwenhuizen |12] and Hartman and Weichman 



fl"3| have considered quantized versions of the spherical model. In both cases, after some 
manipulation, [0,|13|] it is possible to write the Hamiltonian in the form of Eqs. ([L8] - 
|1 . 1 1| ) . However, the constraint considered by Hartman and Weichman corresponds to the 
case a s = 0, while that of Nieuwenhuizen corresponds to a p = 0. As we shall see, this 
difference has important physical consequences. In particular, to exhibit the two conjugate 
ordered phases, it is necessary that both a's be non-zero. 

Effects of quenched disorder can be included by adding a random field to the Hamiltonian 

H dis = J2 h A- (1-13) 

j 

For our purposes, the ensemble of random fields is adequately specified by its second moment, 

ikk>} = k,-k>f(k), (1-14) 

where the [ ] signifies configuration averaging. (We imagine that [hj] = 0, although this 
is unimportant.) Since the random fields are, presumably, generated by remote quenched 
charges, we will typically suppose that f{k) vanishes as k — > as f(k) ~ k 2 , but is otherwise 
positive at all other values of k. (Note that the two previously cited studies |l2| , |i"3| considered 



quantum spherical models with random exchange interactions so as to construct a spherical 
version of a spin-glass; this problem is interesting but considerably more complex than the 
one considered here.) 

We now return to the additional terms (signified by ...) in Eq. ( |1.3|) . The Hamiltonian, 
as written, depends on k only through A(fc). It turns out that this implies a degeneracy 
(that is non-generic for lattice systems) which is a consequence of a lattice version of the 
rotational symmetry of free space and hence would be exact in a continuum version of the 
model. We will lift this degeneracy through the inclusion of the term 

which in position space corresponds to a second neighbor interaction. This term vanishes 
like fc 4 as k — > 0; notice that, to this order, any perturbation consistent with the symmetry 



2 ^{T. v ^)v(h)\ (1-15) 



of the hypercubic lattice can be written as A' A(k) plus the term in Eq. ( 1.15 ). However, 
any portion of the perturbation which depends only on the Laplacian does not lift the 



degeneracy. In this sense, Eq. ( 1.15 ) is the unique, leading order term that breaks the 
"rotational" symmetry of the continuum. For concreteness, we will always consider the 
model with A positive, unless otherwise stated. 

To summarize, in the rest of the paper we will consider the statistical mechanics of the 
classical model defined by the Hamiltonian in Eq. ( |1.3|) with 
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J(k) = JA(k)+Q A(k) 



-X^2V(k a )V(k b )+K 



;i.i6) 



and with the equation of constraint defined in Eq. (|1.7|) . We will then consider the quantized 
version of the same model defined by the Hamiltonian in Eq. (|1.8|) with H qu defined by Eq. 
( P-.llj ) and the equation of constraint by Eq. ( |1.10| ) with a s = a p = 1. Henceforth, we will 
work in units such that Planck's constant, Boltzmann's constant, the energy scale J, and 
the lattice constant are all set equal to 1. 



II. DISCUSSION OF RESULTS 
A. Classical Model 

For Q = 0, the model described is the usual spherical model. As is well known for this 
case, d — 2 is the lower critical dimension; there is no finite temperature transition in 
d = 2 but the correlation length diverges rapidly, like £o = L(X) exp [2n/T], as T — > 0, 
which is qualitatively similar to the behavior of the 0(n) model for n > 2. In d = 3, 
there is a finite temperature transition at a critical temperature T c = A(X) where A(\) 
is a strongly varying function of A with A(0) ~ 0.79, and the correlation length diverges 
as £o = {4:7r/A) [(T — T c )/T c ]~ u with v — 1 as the temperature approaches T c from above. 
The spherical model is somewhat unsatisfactory for temperatures below T c where there is 
a well defined growth of the order parameter, (Sj) = [{T c — T)/^] 1 ^ 2 , but the (Josephson) 
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correlation length is infinite. (This is one of the problems we hope will be fixed when 1/n 
corrections are included. ||) 

For non-zero Q, ferromagnetic order is completely forbidden. However, modulated order 
of one sort or another is permitted. For small Q, it can be easily seen by minimizing J{k) 
in Eq. ( |1.3[ ) with respect to k that the preferred order occurs for k = q, where V(q) = \J~Q. 
Since q is small in this limit, a first approximation to the physics can be obtained by taking 
a continuum limit, in which only the order k~ 2 and k 2 terms in J are retained; in particular, 
this means that the term proportional to A is neglected since it is smaller by a factor of y/Q 
than the first two terms for k ~ q. However, as is well known in the continuum theory of 
smectic liquid crystals, the transverse fluctuations of such density wave order are sufficient 
to destroy the long-range order which, in the context of the spherical model is equivalent 
to the destruction of any finite temperature transition. (A power law phase, as in the 2c? 
XY model or the 3d continuum theory of smectic liquid crystals is impossible for this simple 
class of models.) In the present context, this is a consequence of the existence of a co- 
dimension one hypersurface of minimizing wave vectors which implies insufficient stiffness 
against fluctuations to permit a finite temperature ordering transition in any dimension. As 
discussed above, this piece of continuum physics is reproduced on the lattice when only the 
first two terms in J are retained (even though they are computed to all orders in k). In 
this case the preferred order occurs at k on the hypersurface A(A;) = ^/Q and therefore so 
long as < y[Q < 4d, no finite temperature ordering transition is possible. (For larger Q, 
even if A = 0, the lattice asserts itself, and ordering can occur.) Thus we see, trivially, that 
in the model with A = 0, there is an avoided critical point for all dimensions d > 2, as in 
Fig. 1, in which T c drops from a finite value for Q = 0, to zero for Q ^ 0. 

We now address the question of what happens to this phase diagram when the symmetry- 
breaking term proportional to A is introduced. In this case, the minima of J{k) occur at 
isolated points in the Brillouin zone, and hence finite-temperature ordering is possible at 
non-zero Q in all dimensions d > 2. For Q < 16 and A positive, there are 2d minima which 
occur at k = ±qe a , for a = 1, d, with V(q) = y/Q, and hence correspond to unidirectional 
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striped phases. (The ordering vectors for all Q and all A are listed in Table 1.) In dimension 
d — 3, as shown in Fig. 1, there remains a finite discontinuity in T c as Q — > 0. Specifically, 
for positive A, T C (Q = 0) = A Q + A X A + C(A 2 ) with A Q = A(0) and A 1 = A'(0) ~ 0.074, 
while 

hm T C (Q) = v^AT c (0) 

with _B ~ 0.14245. For d > 3, finite A results in a continuous, although extremely non- 
analytic behavior of T c at small Q: T C (Q) = T C (0)[1 - B d T c (0)\-^(Q^ d - 3 ^ 4 } where B d is a 
dimension dependent number of order 1. (For what it is worth, for 2 < d < 3, there is a 
finite discontinuity in T c , and indeed, T c is finite for Q — 0, but tends to zero in the Q — > 
limit as T c ~ B^V\(Q) {3 - d)/4 .) 

In Fig. 1, we show the phase diagram for the three-dimensional classical model as a 
fuction of Q and T for fixed, positive A (The solid lines in the figure are the phase bound- 
aries computed for A = 1/4, but the qualitative results are insensitive to the value of A.) 
There appear two dashed lines in the small-Q region of the figure which signify crossover 
temperatures discussed below: 

• 7\ marks the temperature at which the frustration becomes significant. At T > Ti, 
correlations behave essentially like those of the model with Q = 0, while for T < Ti, the 
effect of the frustration is to break the incipient ferromagnetic order into randomly oriented 
"domains" of typical size, f, which at T w T C (Q = 0) has magnitude f « (1/<?) 1/4 , and 
which grows slowly as the temperature is lowered. 

• T 2 marks the temperature at which lattice effects becomes important. For T > T 2 , the 
fluctuations are essentially isotropic, while for T < T 2 , the correlation length begins to 
diverge as the ordering temperature is approached, with the same critical exponent as the 
unfrustrated model, and the correlation functions begin to choose preferred orientations for 
the domains, corresponding to the onset of stripe ordering. One additional pathology of the 
spherical model apparent in this phase diagram (which, we expect, is corrected in order 1/n) 
is that there is no commensurate lock-in whatsoever for small Q, i.e. there is no preference 



for stripe ordering wave vectors q which are commensurate with the underlying lattice. In 
Fig. 3 we show the Fourier transform, G(k), of the spin-spin correlation function, 

G(R) = (S 3 S M ) (2.2) 

at various temperatures for a fixed, small value of Q — 1/2 and A = 1/4. One can clearly see 
the differences in the structure of the correlations in the different regimes of temperature. 

The most dramatic manifestation of avoided critical behavior is the existence of these 
sharp crossover regimes. For instance, we associate the remarkable properties of supercooled 
liquids as they become glassy with the behavior of a uniformly frustrated system in the 
temperature range T ^ T\. The structural correlation functions of certain high temperature 
superconducting materials also exhibit behavior which is strikingly similar to that of the 
various low temperature regimes of the present model. It seems that avoided critical behavior 
is surprisingly robust, at least in d < 3. It occurs naturally in the continuum version of the 
model, and survives lattice effects. Preliminary results from the 1/n expansion || show that 
it persists at finite n. This is a new construct which is likely to find applications in a variety 
of other arenas. 

B. Effects of Quantum Fluctations 

In the quantum model, there are two distinct types of possible ordered states: charge- 
ordered states, in which Sg develops a non-zero expectation value, and superconducting 
states, in which Pg develops a non-zero expectation value. At zero temperature, the ground 
state will be ordered if s exceeds a parameter-dependent critical value s c while the ground 
state is quantum disordered if s < s c . (When factors of % are restored, s 2 — > s 2 /h, in 
Eq. ( |1.10| ) so s 2 is the natural measure of the importance of quantum fluctuations.) The 
superconducting state can occur only for K greater than a critical value, Kq, while the 
charge-ordered state can occur only for K < K ; coexistence, i.e. supersolid order, can 
occur only for the special case of K = K . ( For Q < 16, K = —2\fQ\ more generally, K 
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is that value of K for which the minimum over all k of J~(k) in Eq. ( 1.16| ) is equal to 0.) As 
in the classical model, the transverse quantum fluctuations in the continuum approximation 
(or for A = 0) are sufficient to destroy any possible charge order for small Q, even at zero 
temperature. However, for non-zero A, not only is there a finite value of s c for d > 1, but 
it is a continuous function of Q, even in the limit Q — > 0. Indeed, s c is a rather weak 
function of all parameters, and is always greater than, but approximately equal to 1. (The 
dependence of s c on K and Q for W = 1 is shown in Fig. 4.) Thus, at zero temperature 
and for s > s c , there is a superconducting to charge- density wave transition that occurs 
as a function of K. Similarly, as a function of decreasing s, there is a superconducting to 
quantum-disordered transition which occurs for K > Kq, and a similar charge- density wave 
to quantum-disordered phase transition for K < K . 

At finite temperatures, the behavior of the system for K < K and s > s c is similar to that 
of the classical frustrated model described above, while for K > K , the superconducting 
to normal transition occurs in qualitatively the same way as in the ordinary ferromagnetic 
spherical model. However, so long as \K—K \ is small, both the normal and superconducting 
states will exhibit substantial correlations which resemble the nearby charge-density wave 
phase. The same two crossover temperatures that appear in the correlation function of the 
classical model calculations (Fig. 3), but in contrast to that case, the spin correlation length 
saturates at a long but finite value for T < T c where T c is the superconducting transition 
temperature. 



C. Effects of Disorder 

For Q > 0, disorder destroys the possiblity of a charge-ordered state in all dimensions 
d < 4. For K < K , the ground-state in d < 4 has more or less extended short-ranged 
charge correlations, depending the strength of the disorder, but no true long-range order. 
(Indeed, the behavior of the disordered system in dimension d is qualitatively similar to 
that of the ordered system in dimension d — 2, or in other words the disorder produces the 
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standard "dimensional reduction" associated with random field problems fL7j ) For Q = 0, 
however, the fact that the mean-squared random field, f{k), vanishes at small k insures that 
a finite transition temperature to a charge-ordered (ferromagnetic) state survives up to a 
critical magnitude of the disorder. Thus, in all dimensions d < 4, the model has an avoided 
critical point. By contrast, weak disorder has relatively little effect on the ground state 
and low-temperature properties of the superconducting state unless K — K is quite small. 
For fixed disorder, and K > K Q , the superconducting T c tends continuously to zero as K 
is decreased, and always vanishes before K > Kq. (Recall that, in the absence of disorder, 
the superconducting state would have given way to a charge-ordered state, both with finite 
transition temperatures, at precisely K = K .) 

In Fig. 5, we show the phase diagram as a function of temperature and position along the 
representative trajectory through parameter space indicated by the arrow in Fig. 4, both in 
the presence and absence of disorder. 

D. Some Exact Ground-states for the Ising Model 

An amusing side benefit of the present analysis is that, for certain ranges of parameters, 
the exact ground states of the classical spherical model are also the exact ground states for 
all n of the 0(n) model with the same Hamiltonian, including the n = 1 Ising case. (These 
are listed in Table I.) The proof is simple: For the spherical model, the ground states are 
found by minimizing the Hamiltonian with respect to spin configurations, subject to one 
global constraint. For the Ising model, the same minimization problem must be solved, but 
now subject to iV — 1 additional constraints to insure that each spin has length 1. Thus, the 
ground-state energy of the Ising model must always be greater than or equal to that of the 
spherical model. If it so happens that the ground state of the spherical model is an Ising 
state, then this state must be a ground state of the corresponding Ising model. To extend 
the proof to the 0(n) model, we simply consider n copies of the spherical model with one 
coupled constraint. 
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In this way, we can prove that there exists a region of parameter space for which the 
corresponding frustrated 0(n) model has the spherical model ground state (and is thus 
independent of n). In three dimensions these states are: 1) a six- fold degenerate width 1 
stripe with q = (7r, 0,0), 2) a six-fold degenerate columnar state with q = (7r, 7T, 0), and 
3) a two-fold degenerate Neel state with q = (11,71,71). (The degeneracies apply to the 
Ising model; for n > 1, the 0(n) models have an additional continuous degeneracy.) In 
addition, it is possible to prove that there exists at least a surface in parameter space on 
which the ground-state is: 1) a 12-fold degenerate width 2 stripe with q = (tt/2,0,0), 2) a 
24-fold degenerate rectangular columnar state with q = (tt, tt/2, 0), 3) a 12-fold degenerate 
rectangular Neel state with q = (71,71,71/2), and 4) a 16-fold degenerate width 2 Neel state 
with q = (ti/2,7i/2,7i/2). (Surely, for the Ising model there is a commensurability energy 
which stabilizes these latter states in a finite region of parameter space.) Similar results can 
straightforwardly be obtained in other dimensions. 



III. METHOD OF SOLUTION 



A. The Classical Model 



Consider a model of the form of Equation (|1.3| ) where for convenience, we add a counter- 
term as in Eq. ( |1.12| ) defined so that the minimum value of J is zero. The standard version 
I H| of the spherical model dictates that we integrate over all configurations (S^) subject 



to the constraint Eq. ( |1.7|) . A simpler solution, that was introduced as early as 1952 Jt|, 
employs the method of Lagrange multipliers: The original Hamiltonian is augmented by the 
term ^J2k \^k\ 2 and the integration takes place over all finite energy spin configurations. 
The model is now unconstrained and quadratic, so all quantities can be computed readily. 
As long as /i > 0, this can be done without apologies and equation of constraint becomes 
an implicit equation for fi(T): 

f = ^ )s ^mr, (3 - 1) 
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or, in the thermodynamic limit, 



1 



d d k 1 



(3.2) 



T 



and the integral is over the first Brillouin zone. If this equation cannot be satisfied for any 
value of fi, we are at or below criticality. Since $ is a monotonically decreasing function of 
T c is determined according to 



For T < T c , the Lagrange multiplier is set to zero. The total population of the finite modes 
therefore is deficient, and the remainder is identified as a condensate to be distributed among 
the zero modes. 

It is intuitively clear that at least as far as the equilibrium properties of local observables 
are concerned, the above procedure is equivalent, in the thermodynamic limit, to the original 
constrained model. Such results have been established with a large degree of generality - 
more than sufficient to cover the cases of interest here. See, for example, [Tj| (in particular 
Theorem 1) and references therein. 

Most of the claims made in the earlier sections thus amount to explicit calculations or 
elementary analysis of the function J{k). In particular, the internal energy per site U as a 
function of temperature is simply 



and all other thermodynamic quantities can be determined by taking appropriate partial 
derivatives of this expression. For T < T c , the magnitude, m, of the condensate is simply 



The spin-spin correlation functions can be straightforwardly calculated according to 




(3.3) 



U = \ [T - M (T)] 



(3.4) 



m 2 = 1 - T$(0) = [T c - T}/T c . 



(3.5) 



G(k,T)^(\S k \ 2 )=T/[J(k) +f i(T)}. 



(3.6) 
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as is easily seen using the so called method of generating functions: Add small terms to 
the Hamiltonian (well localized in position or momentum space) and differentiate the free 
energy with respect to the appropriate coupling. From Equation ( |3.6| ), the inverse of the 
correlation length may be extracted: 

UT)- 1 = mm\lm{q° a }\ (3.7) 

(in general, may be anisotropic, hence it is labeled by a direction) where q° is the 
solution of the implicit equation 

J(q°)+fi(T) = 0, (3.8) 

and the minimization in Eq. ( |3.7|) is over the set of solutions with all other components of 
q real. (In general, may be anisotropic, hence it is labeled by a direction.) 

Crossover temperatures can, of course, never be determined from a single sharp criterion. 
However, near the avoided critical point, a crossover occurs in a narrow range of tempera- 
tures, so it is both useful and appropriate to define a crossover temperature explicitly. The 
crossover temperature, Ti (below which the frustration becomes "important") is therefore 
defined as the solution of the implicit equation 



K T i) = 2 \/Q- ( 3 - 9 ) 

where the factor of 2 is chosen for aesthetic reasons. Notice that for an avoided critical 
point, Ti — > T C (Q = 0) as Q — > 0. Along a trajectory in the Q — T plane which lies at fixed 
distance above T\(Q), the correlation length approaches a finite limit as Q —>■ 0, while below 
Ti, the correlation length diverges in this limit. The lower crossover temperature T 2 (below 
which the lattice anisotropy becomes "important") will be defined - with the same degree 
of arbitrariness as T\ - as the solution of 

H(T 2 ) = XQ. (3.10) 

For small Q and A and Ti 3> T ^> the structure factor is dominated by a sharp ridge at 
| A; | = q, while for T 2 > T, there are six sharp peaks at \k\ = q and k along a coordinate axis. 
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With our definition, at T = T 2 , the modulation in the magnitude of S(k) (as a function of 
angle) with fixed \k\ = q, is comparable to the height of the ridge. 

B. The Quantum Model 

The quantum model is solved in much the same way as the classical model, through the 
introduction of a chemical potential to enforce the constraint, which reduces the problem to 
a set of decoupled harmonic oscillators. Thus, for all temperatures above T c , /i is implicitly 
determined from the relation 



s 2 = $ gn (/i,T) 



(3.11) 



where 



% u (fi,T) = - 



d d k 
(2vr) a 



WA(k)+fj, | J(k)+fj, 



2n{ut/T) + 1 



(3.12) 



WA(k)+/i J(k)+/i 



(3.13) 



and n(x) = [e x — 1] _1 is the Bose occupation factor. (The behavior of the dispersion relation 
as k — >• is somewhat peculiar, but this has no effect on the present results. [|14|| ) Again, $ 
is a monotonically decreasing function of fi so, at T = 0, the critical value of s is determined 
from the equation 



s c 2 = $(0,0) 



(3.14) 



and, when s > s c , the critical temperature T c is determined according to 



s 2 = $(0,T c ). 



(3.15) 



For s > s c and T < T c , the magnitude of the condensate, m, is determined according to 



m 2 = S 2 -$(0,T), 



(3.16) 
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where, as discussed above, m is the magnitude of the superconducting condensate for K > 
K Q and the amplitude of the charge- density wave for K < K . The dynamic spin correlation 
function (since in quantum mechanics, dynamics and thermodynamics are intimately related, 
the Fourier transform of the two-time spin correlation function is the fundamental quantity) 
is easily seen to be 



G{k,u) 



WA(k)+/j 



2 \ J(k)+fi 
while the static structure factor, 



{ [n(w s /T) + l] 6(u - u % ) + n{^/T)5{u + , (3.17) 



G(k) = I —G(k,uj) = —a = 



WA(k)+/2 r 



2n{ut/T) + 1 



(3.18) 



The inverse charge-ordering correlation length (which, again, will generally be anisotropic) 
is 



(T) = min |lm{gf } 
where <f h is the solution of the implicit equation 

J{H») + H(T) = 0, 



(3.19) 



(3.20) 



and the minimization is again performed with all other components of k real. A similar 
expression for the superconducting correlation function is easily obtained by inverting the 
term in the square root in Eq. ( |3.17|) and a superconducting correlation length obtained 
from the solution of the implicit equation 



WA{f c ) +/i(T) = 0. 



Finally, the internal energy per site U as a function of temperature is simply 



U = U (T) - -n(T) 



where 



Uo(T) 



d d k 



2n{ut/T) + 1 



(3.21) 



(3.22) 



(3.23) 



(2yr) d k 

is the internal energy of independent harmonic oscillators, and all other thermodynamic 
quantities can be determined by taking appropriate partial derivatives of this expression. 
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C. The Model with Disorder 



Because of the harmonic nature of the model, the effect of an arbitrary configuration of 
random fields {hj} can be formally accounted for by shifting all spins according to 



(3.24) 



for either the classical or the quantum model. The result is an additive term to the internal 
energy, 



d d k 


h 


2 


(27r) d 





U -> - 



and a shift in the implicit equation for fi brought about by the substitution 

2\ 



(3.25) 



Si 



s, 



hi 



2 ' 



Upon configuration averaging, this relation implies 

d d k f(k) 



(3.26) 



® dis (li,T) = $ gu (/i,T) + J 



(2vr) c 



2 ! 



(3.27) 



(and the obvious corresponding equation for the classical model). The spin correlation 
function in the presence of disorder is altered both by the implicit change produced by 
the altered temperature dependence of /i, and by the addion of a zero fequency additive 
contribution 



G dis (k,u) = G(k,u) + 



*M/(fc) 



2 ' 



(3.28) 



J{k)+n 

It is clear from Eq. ( |3.27|) that for K — K negative & dis diverges as \i — > — K for all 
dimensions d < 4, so that no charge ordered state is possible. Moreover, UK — K is 
positive but small, $dj S (0, T) will have a large, additive contribution from the disorder term, 
from which it follows that the superconducting T c must always vanish as K — > K + AK 
where AK is a positive, increasing function of the strength of the disorder. 
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IV. MATHEMATICAL DUMP 



We now provide the calculational details which underly the various claims made in the 
previous sections. Our starting point is an elementary calculation concerning the stationary 
points of J(k): 
Proposition IV. 1. Let 

J Q (k) = E V ( k «) + V(k a )V(k b ) + g + K (4.1) 

— * 

with V{k a ) = 2(1 — cosk a ) and A ^ 0. Let p be a stationary point of J{k) in the first 
Brillouin zone, — tt < k a < tt. Then any component of p may be or tt but the components 
that are not or tt are equal in magnitude to each other. In particular, if p is a stationary 
point then \p a \ = 0, tt, or q^ )S with g^ jS satisfying 

l+m+( s-m^)\= {u J VMf (4.2) 

where I is the number of components equal to tt and s is the number of the remaining 
components that are non-zero. This theorem applies in arbitary dimension, d > 1. 
Proof. Upon differentiating J{k) we find 

[1 + A E V(k b ) - Q ] sin k a = 0. (4.3) 



This is obviously satisfied if k a = or tt. Now suppose that k satisfies Eq. (|4.3| ) with two 
or more components, k b and k c not or tt. Then, subtracting, we get 

XV(k b ) = XV(k c ) (4.4) 

which implies \ki,\ — \k c \. If there are I directions where |/c a | = tt and s directions where 
\k a \ is not or tt then the magnitude of these remaining components satisfies the stated 
equation. □ 

What follows is the starting point for both the analysis of the low temperature behavior 
in the system as Q — > and for the analysis of the ground state space. For the most part, 
the indices on the q defined above will be understood from context and omitted. 



Proposition IV. 2 Let J(k) be as described in Proposition IV. 1. Then for A > ; and Q < 

16, the minimizing wave vector has a single non-zero component of magnitude q satisfying 



V(q) = sjQ. (4.5) 

Proof. If A = 0, the function has an absolute minimum of 2y/Q which is achieved if 
J2 a V{k, a ) = \[Q. Provided that \[Q < 4, this value can be obtained even if A > by 
a vector with only a single component that has the above stated magnitude. □ 

We are now ready for our principal result for this section: 
Theorem IV. 3. In three dimensions, for the model described in Proposition IV. 1, there is 
an avoided critical point at Q = 0. In particular, for fixed A > ; let T c (0) (= T c (0; A)) be 
the critical temperature for the (Q — 0) ; ferromagnetic version of the model: 

1 1 r d 3 k 

Tj0) ~j2^J\ ka \<, E a V(k a ) + \XE a ^V(k a )V(k b ) (4 ' 6) 

and T C (Q) given as in Eqs. (P-Q) and p.^) with Jq{k) given in Eq. (fJTTL) with K = K = 



-2^Q. Then 



< hm = — — + -=, (4.7) 



T c (0) Q-oT c (Q) T c (0) 

Where B = f dQ[sm 4 6 cos 2 0sin 2 + cos 2 9 sin 2 6]~? w 0.14245 
Proof. For Q 1, let q ~ Q 3 denote the solution of 



V(q) = (4.8) 

and let A be a small number independent of Q, the precise specifications of which will be 
detailed later. (In most of what is to follow, the distinction between q and is practically 
irrelevant - q may simply be regarded as convenient notation for Q^.) It is clear that for the 
large-fc portion of the integral, \k\ > A, as Q — > 0, nothing particularly interesting happens 
to the integrand [Jq 1 ^)}. Hence 

f d 3 k f d 3 k 
lim / — = / —. 4.9 

Q-0i|fe|> A J Q (k) J\k\>A J Q=0 (k) 
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Notice that, for small A, the right hand side is only just shy of 1/T C (0). 

Now, let < t < 1 and define the two constants A± = 1 ±g*, which have the property that 
A± — > 1 as Q — > 0. In terms of these, we break the remaining region of integration in three: 
1) A small k region, TZi with < \k\ < qA^; 2) A critical region, 712 with qA_ < \k\ < qA + ; 
3) An intermediate k region TZ 3 with qA + < \k\ < A. We shall show that, as Q — > 0, the 
contribution from region 1 vanishes, the contibution from region 3 plus the contribution 
from large k in Eq. ( |4.9|) converge to 1/T C (0), and the contribution from the critical region 
approaches B / VX 

• In 1Z\, as an upper bound, we will neglect the A-perturbation in the denominator. Noting 
that in the specified region, J2 a V(k a ) + ^ ® . . is increasing and that, in general, V(k a ) < 
kl, we have 



<47rg 4 / ^-3—. (4.10) 

J\k\<aA. (k - OlV 



iKi J Q (k) Hi k 2 + Q/k 2 -2y/Q J\k\< q A^ (k - Qr 

The last term is bounded by constants times g( 3 ~*) as Q — » and thus the limiting contri- 
bution from this region may be neglected. 

• Next, consider the intermediate region IZ3. Since Jq = JTo + Q/ J2 a V{k a ) — ^\fQ it follows 
that 

d 3 k _ f d*k + 1 [2y/Q-Q/E a V(k a )]d 3 k 



lK3j Q (k) Hj (k) H J (k)J Q (k) 
Now the first term on the right hand side of Eq. ( [4.1 1|) may be added directly to the right 



hand side of Eq. (|4~9 ) to obtain, in the small Q limit, exactly 1/T C (Q = 0). Let us show 
that the second term vanishes as Q — > 0. In this region, the numerator is positive so we may 
discard the term involving Z) V(A; ). The denominator is made smaller if we set A = so 
let us do that as well. Let D denote any positive constant for which V(k a ) > Dk 2 holds for 
all k a with \k a \ < ir. Putting these together, and canceling powers of k the task is to show 



that 



Q h Y.aV{k a )-2^Q + Y; ^- ) (4 ' 12) 
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vanishes in the Q 



limit. For the term 



in the denominator, we may use the 



bound V{k a ) < k\ and, further, for the other appearance of J2 a V(k a ), we may replace 
this sum with k 2 — Ek 4 where E is some constant independent of Q (or A). With these 
estimates in tow, the problem is one dimensional and, were it not for the quartic term, would 
be entirely trivial. In any case, we are now reduced to showing that 



vanishes as Q — > 

Let us now assert that A has been chosen small enough so that (for all q sufficiently small) 
throughout the range qA + < k < A, 



Indeed, provided that q is sufficiently small, the inequality clearly holds at the lower limit 
and, upon comparison of derivatives, the desired inequality holds throughout the entire 
range provided that EA is somewhat less than one. What is left after these estimates can, 
essentially, be done by hand: 



which indeed vanishes as Q — > 0. 

• Thus we are left with the critical region 7Z 2 - In this region, the deviation of k from q 
is so slight that we are essentially in the position where we can "expand and neglect". In 
particular, upper and lower bounds may be derived, in this region, by setting various items 
to their maximum or minimum value. Since the procedure is similar on both sides, we will 
be content with an explicit derivation of a lower bound on the remaining integral that agrees 
with the stated formula. As will become apparent, an upper bound can be derived in the 
same fashion. 

In what follows, all constants C n will be functions of Q with the property that C n — > 1 
as Q — > 0. It is slightly easier to work with the quantities V(k a ) instead of the k a hence we 
define the variables 




(4.13) 




(4.14) 




(4.15) 
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v a = iV(k a )]z (4.16) 

and note that dv a = 1 2s in'(lfc ) l — Cidv a . It is further clear that in the region \k — q\ < q 1+t , 
we may write d 3 k > C2q 2 dvdQ. Next, it is noted that the image of the region \k — q\ < q 1+t , 
in f-space, also contains a spherical shell of size C^ 1+t ^ . We will confine our attentions to 
this smaller region. Thus, 

d 3 k ^ „ r dudVL 



> C 2 ( <iud\l 



where u = v — q, 



R(Q) = lim \ ^Yl v l v l = N n2 9 cos2 9 + sin4 9 cos2 ^ sin2 $ ■ ( 4 - 18 ) 

For R(il) 7^ 0, the u integral is easily preformed to yield 

d 3 k > Cfe r (4 lQ) 



n 2 J(k) y/\J jR(^j 



where 



a Q = tan 1 



2C 3 q t - 1 /^/\R(n) 



tt/2. (4.20) 



Hence, it is seen that 



f d 3 k n f dfl , 

lim / > — -= / . 4.21 

Q-oJn 2 J(k) 2^\J ^fR(a) 

as claimed. 

An upper bound follows almost the identical derivation with a renaming of the constants 
C n . □ 

Corollary. In dimension d > 3, there is a "nearly avoided critical point" for the stated 
model in the sense that for fixed A > 0, as Q — > 0, 

' '' ~Q*r\-*B d (4.22) 



T C (Q) T c (0) 

where B d = (27r)~ d j* n dtt[R d (tt)]-^ with R d (Q) = \Y,a^bX a x b U*=i- 
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Proof. If we use the same division of regions and follow step for step the analysis of Theorem 
IV. 3, we see that all the estimated quantities that go to zero are now multiplied by an extra 
factor of fc( d_3 ) (from the volume element). This always amounts to an extra factor of g( d ~ 3 ) 
Similarly, the term that was of primary interest works in a similar way but is multiplied by 
q(d-3)^ Everything else that was finite produces ^4— mutatis-mutandis. In particular, we 
get 

feeder ^k 1=riBj (423) 

□ 

Our discussion of the ground state space now picks up where Proposition IV. 2 left off. 
Right now, for A > and Q < 16 the situation is well under control and these are the sorts 
of results that we seek throughout the phase plane. For the sake of brevity, we will focus 
our attention on the cases of principal interest, namely d = 3 and A > 0. 

Proposition IV. 1 tells us that in any region, there are only a finite number of possibilities 
to consider; however, for d = 3, this turns out to be eight additional distinct modes (other 
than (q, 0,0)) and five new regions. Notwithstanding, we will attempt to be as brief as 
possible and still lay claim to a rigorous proof; this is most efficiently carried out by writing 
out five separate sub-propositions. A complete list of the competing modes as well as the 
regions of interest can be found in Fig. 2. In the up and coming, the various state (or modes) 
will not be distinguished from their reflection or coordinate axis exchange equivalents. Thus, 
e.g. Proposition IV.4.i below really pertains to six ground (equivalent) states. 
Proposition IV. 4. i. In the region A > 0, Q > 16 and Q < 16(1 + 4A) ; the ground state is 
of the form (tt, 0, 0). 

Proof. Since Q > 16, the (q, 0, 0) states need not be considered because there is no solution 
to Eq. ( |4.2| ) for q. Similarly, for the mode (it, q, 0), the defining equation reads (1 + 4V^) 2 = 
(where here, and in what is to follow, we use the notation V q = V(q)). Hence this 
mode is forbidden if Q < 16(1 + 4A). In the same region, for the same reason, the mode 
(tt, q, q) is disallowed and similarly, the mode (tt, tt, q) is forbidden for Q < 64(1 + 8A). Now 
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the "energy" for the current state, <7(ir, 0, 0), is simply 4 + Q/4. In the stated region, this 
is a whole lot less than 8 + Q/8 + 16A and 12 + Q/12 + 48 A which eliminates (tt, it, 0) and 
(7r, 7r, 7r) from consideration. This leaves as contenders only the modes (q, q, q) and (q, q, 0). 
The energy for (q, q, 0) is given by 

J { q, q, 0) ee £ M)0 = 2V, + ^ + AV g 2 . (4.24) 

Now in general, the energy of any state increases with A however £^0,0 is independent of A. 
It is therefore sufficient to establish £ q ^o > £n,o,o at the lower-right boundary of the region 
under consideration. To this end, we write 

S q , q>0 = 2V q + ^- + 4X[2V q ] - 16A + A(4 - V q f. (4.25) 

L V q 

Neglecting the quadratic term and minimizing at Q — 16(1 + 4A), this is maximized when 
2V q = 4 and weighs in at exactly 4 + Q/4. Using the fact that we are interested in Q's that 
are strictly less than 16(1 + 4A) and, the fact that £ qA ^ is a strictly increasing function of 
A we find £ q , q ,o > £n,o,o- The mode (q,q,q) has energy 3V q + Q/3V q + 3AV g 2 which (when 
regarded as a function of 3V^) is manifestly larger than £ qA $- 

Proposition IV.4.ii. In the region A > ; 16(1 + 4A) < Q < 64(1 + 4A) ; the ground state 
is of the form (n, q, 0) . 

Proof. Notice that this is exactly the region where Eq. ( f4.2|) has a solution for the mode 
(7r, q, 0). The mode (71", it, q) has no solution in this region (and nor does (q, 0, 0)). We write 

£ nq0 = mm\uj + — + 4Acj-16A1 (4.26) 



(with 00 ee 4 + V q ; note that the minimizing 10 equals JQ/(1 + 4A)). On comparison to £ q , q ,o 
as expressed in Eq. ([4.25|) , it is clear that £ n , q ,o is lower. This also eliminates (q, q, q) on the 
basis of the final argument in Proposition III.4.i. Next we have £ nt7Tt o = 8 + Q/8 + 16A and 
plugging uj = 8 into the expression on the right-hand side of Eq. ( f4.26| ), this is exactly what 
we get. The minimizing u will do better. 

Thus we may turn our attention to the mode at (it, q, q) - the ones at (tt, 0, 0) and (n, it, 7r) 
then follow immediately. The energy, £-K^ q admits the expression 
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= 4 + 2V q + + 4A(4 + 2V q ) - 16A + 8XV q 2 . (4.27) 



Regarding this as a function of 4 + 2V q , and comparing to Eq. (|4.26| ), this energy is clearly 



larger than £^,0 whenever we are in a region where a minimizing u for Eq. Q4.26 ) exists. 



Proposition IV.4.iii. In the region A > 0, (1 + 4A) < Q/6A < 64(1 + 8A) ; the ground 
state is of the form (it, it, 0). 

Proof. The modes (q, 0, 0) and (71*, q, 0) are eliminated from consideration. Similarly, V q for 
(q, q, 0) must satisfy 

Q 
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but 



Q > — 0— > 64 (4.29) 



1 + XV q ~ 1 + 4A 
so for (q, q, 0), cannot get big enough. 

For the mode at hand, the energy £^,^,0 is given by the simple formula 

8*^,0 = 8 + Q/8 + 16A. (4.30) 

Subtracting this from various expressions for the energy of various other modes, we see 
8 7r fl,o ~ £w,ir,o oc Q — 16(1 + 4A) > 0, S^ n>1t — = 4(1 + 8A) — Q/24 > and similarly, 
£-7T,n,q — £tt,tt,o = V q (l + 8A — Q/(8)(8 + V^)) > 0. Notwithstanding the losing status of the 
state (it, 7t, q), for the benefit of the final two states on the list, let us write its energy: 

8 m = 8 + V q + -0— + A[16 + 8V q ] = r) + 5 + 8X V - 48A (4.31) 

0+ V q 1] 

where 77 = 8 + V q . Similarly, £ q , q , q may be written 

£ - a + 9. + 8Acr - 48A + A[-a 2 + 48 - 8a] (4.32) 
a 3 



with a = 3V q . (Of course the g's referred to in Eqs. ( (4.31| ) and ( |4.32| ) pertain to different 



solutions of Eq. (|4.2|) and are not to be identified with one another.) It is easily checked 
that !<r 2 + 48 — 8a > for a < 12 and hence it is seen that £ q , q , q > £-k,-k, q - 

25 



Similarly, we write 



^ = 4 + 2% + + X(8V q + V 2 ) = ( + ^ + 8AC - 48A + \[-( 2 - 6( + 36] (4.33) 

and again the term in square brackets is positive if ( < 12. 

Proposition IV.4.iv. In the region A > and 64 < Q/(l + 8A) < 144, the ground state 
is of the form (ir, ir, q). 

Proof. The modes (q, 0, 0), (ir, q, 0) and (q, q, 0) are already out of the loop. In Proposition 
IV.4.iii, we have (just) shown £ q , q , q > £-K,-K, q and £n, q , q > £ir,n, q so these are out. Now at the 
minimum, V q satisfies (8 + V q ) 2 = Q/(l + 8A). In the specified region, this is not solved by 
V q = or V q = 7T so evidently > £-w,-K,q and £^,0 > £-K,-w, q - Finally, we dispense with 

(ir, 0,0) by noting that in this region, £^,0 < £-k,o,o- 

Proposition IV. 4. v. In the region X > 0, Q > 144(1 + 8A) the ground state is of the form 

(7T,7T,7r). 

Proof. Examining Eq. ( |4.2j ) for V qi in the state (gi, q^, 93), we have 

1 + ^ + ^ = (K tl+ ^ + v.). - (434) 

The right hand side is larger than Q/144 and the left hand side is smaller than (1 + 8A). 
Thus there are no solutions with any component not equal to or ir. It follows from earlier 
results (or it can be easily checked) that in this region, £^ <K{K < £^,0 < £n,o,o- 

□ Corollary. Along the line = 36, the ground state is (ir, |,0) and along the line 
jQ-^ = 100, the ground state is (ir,ir, |) 

Proof. This follows from setting V q = 2 in the appropriate region and solving Eq. (f4.2|) for 

Q(A) 
□ 

Remark The case A < is far easier to analyze. Indeed, writing J(k) = A(k) + ||A|[A(A;)] + 
Q/A(k) — ||A] J2a V(k a ) 2 , let k denote any wave vector, and let q satisfy dV(q) = A(k). It is 
clear that J(q, . . . , q) < J(k) unless k = ((q, . . . , q)). Indeed, this amounts to showing that 
J2 a V 2 (k a ) > i[A(A;)] 2 and the latter is just (the discrete form of) Holder's inequality which 
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holds as an equality if and only i f the V(k a ) are independent of a. Evidently the minimizer 
is "diagonal" fr om which it is easy to see that (in three dimensions) the Neel ground state 
dominates in the region Q < 144 — 1152|A| for A < 0. 
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A > 


< Q < 16 


q=<q,0,0> 


V(q) = 

\ 1 / v ^/ 


A > 


16 < Q < 16(1 + 4A 


q =< vr,0,0 > 




A > 


16(1 + 4A) < Q < 64(1 + 4A) 


q =< 7r,g,0 > 


F(g) = v /Q/(l + 4A)-4 


A > 


64(1 + 4A) <Q< 64(1 + 8A) 


q =< ir, 71, > 




A > 


64(1 + 8A) < Q < 144(1 + 8A) 


q =< vr,7r,g > 


V(q) = yjQ/(l + BX)-8 


A > 


144(1 + 8A) < Q 


q =< 71,71,71 > 




-1/8 < A < 


< Q < 144- 1152 |A| 


q=< q,q,q> 


complicated 


A < 


144 — 1152 |A| < Q 


q =< 71,71,71 > 





Table 1: Values of the three dimensional wave vector, q, that minimize J{k) in different 
ranges of Q. This determines the locations of the peaks in the structure factor in the 
disordered phase and the location of the Bragg peaks in the charge-ordered phase. The 
value of q labeled "complicated" is given by the solution of V(q) = v where v is the solution 
of the cubic equation 9v 2 = Q + 18|A|v 3 . 



VI. FIGURES 

1) Phase diagram in the T - Q plane for the classical model in d = 3. 

2) Pictures of the various ordered phases of the Ising version of the model for positive A. 

3) G(k) for Q=l/2 and A = 1/4 
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4) Contour plot of s c for fixed A = 1/4 and W — 1. The line with the arrow is the 
trajectory through parmeter space referred to in Fig. 5. 

5) The phase diagram as a function of temperature and position along the representative 
trajectory through parameter space indicated in Fig. 4, both in the presence and absence 
of disorder. 
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